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ABSTRACT 

Nuclear reactions occurring at densities ~ 10 12 g cm -3 in the crust of a transiently accreting 
neutron star efficiently maintain the core at a temperature ~ (5-10) x 10 7 K. When accretion 
halts, the envelope relaxes to a thermal equilibrium set by the flux from the hot core, as if 
the neutron star were newly born. For the time-averaged accretion rates (<^ lO~ lo M yr -1 ) 
typical of low-mass X-ray transients, standard neutrino cooling is unimportant and the core 
thermally reradiates the deposited heat. The resulting luminosity is ~ 5 x 10 32 -5 x 10 33 ergss -1 
and agrees with many observations of transient neutron stars in quiescence. Confirmation of 
this mechanism would strongly constrain rapid neutrino cooling mechanisms for neutron stars 
(e.g., a pion condensate). Thermal emission had previously been dismissed as a predominant 
source of quiescent emission since blackbody spectral fits implied an emitting area much 
smaller than a neutron star's surface. However, as with thermal emission from radio pulsars, 
fits with realistic emergent spectra will imply a substantially larger emitting area. Other 
emission mechanisms, such as accretion or a pulsar shock, can also operate in quiescence and 
generate intensity and spectral variations over short timescales. Indeed, quiescent accretion 
may produce gravitationally redshifted metal photoionization edges in the quiescent spectra 
(detectable with AXAF and XMM). We discuss past observations of Aql X-l and note that 
the low luminosity (< 10 34 ergss -1 ) X-ray sources in globular clusters and the Be star/X-ray 
transients are excellent candidates for future study. 

Subject headings: accretion, accretion disks — stars: neutron — stars: individual (Aquila X-l, 
SAX J1808.4-3658) 
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Introduction 



blackbody fits (as in radio pulsars, Rajagopal fc Roman: 



1996t ^avlin, Pavlov, fc Shibanov 1996| ). This alleviates 



The temperature of an accreting neutron star's core 
and crust is a subject of much current interest, both 
for magn etic field evolution stu dies and thermonuclear 
burning ( Brown fc Bildsten 1998 and references therein). 
For neutron stars (NSs) that steadily accrete at M ~ 
10 _11 -10 -9 Mq yr" 1 , balancing the recurrent heating from 
thermally unstable hydrogen/helium burning with the 
cooling from neutrino emission and radiative dif fusion re- 
quires cor e temperatures T c — (1-2) x 10 s K (Ayasli & 
Joss 1982; Hanawa & Fujimoto 1984). For these accre- 
tion rates, no more than a few percent of the hydro- 
gen/heliu m burning luminosity diffuses inwards and heats 
the core (Fujimoto et al. 1984, 1987). There is an ad- 
ditional furnace in the inner crust (densities above neu- 
tron drip, ~ 5 x 10 n gcm~ 3 ), where the compression 
of matter by accretion induces electron captures, neu- 
tron emissions, and pycnonuclear re actions that release 

a 



lMeV/mp 10 18 ergsg _1 ( [Hacnscl fe Zdunik 



1990). These reactions send most of their heat into the 
core and therefore play a more important role in maintain- 
ing the interi or thermal balance of a constantly accreting 
neutron star (Miralda-Escudc, Paczynski, fc Hacnscl 1990 ; 



Bi ldsten fc Brown 1997| ; prown fc Bildsten 1998^ 

For the case of transient accretors, the heating of the in- 
terior is not so simple. Most of the heat released by hydro- 
gen/helium burning in the upper atm osphere leaves imme- 
diate l y during the unstable burning ( Hanawa fc Fujimotcj 



1986; Fujimoto et al. 1987). Moreover, the time between 
accretion outbursts is much longer than the cooling time of 
the NS atmosphere. As a result, H/He burning alone can- 
not heat the core to the interior temperatures of a steadily 
accreting star. This is an important point, since thermal 
emission from transiently accreting NSs would be observ- 
able when accretion halts , if the core were sufficiently hot 
( van Paradijs et al. 1987; |Verbunt et al. 1994 ; Asai ct al. 
1996b| ; ICainpana et al. 1998a| ). 

In this Letter, we show that direct nuclear heating 
in the deep neutron star crust heats the interior much 
more strongly and naturally explains a large fraction of 
the quiescent emission seen from transient NSs. During 
each outburst, the heat released from the crustal reac- 
tions flows into the core and maintains it at T c « 10 s K 
(Bildsten & Brown 1997). Once the core is in steady- 



state, the heat radiated during quiescence must equal 
the fraction of the heat deposited during the outburst. 
As we show in § 2, in quiescence the NS then emits 
a time-averaged luminosity ~ (1 MeV/m p )(M) ~ 6 x 
10 32 ergss _1 ((M)/10 -11 M Q yr -1 ). For time-averaged (we 
here mean over the recurrence interval) accretion rates of 
(M) <J 10~ 10 M Q yr -1 , this amount of thermal emission 
is unavoidable unless the neutrino cooling is much more 
rapid than the standard mechanisms we considered. 

Of the NS transients detected in quiescence, most have 
a luminosity L q ~ 10 33 ergss" 1 , and the ratio of outburst 
to quiescent fluence is constan t to within a factor of ~ 10 
(Chen, Shrader, & Livio 1997). The rough agreement be- 
tween our estimate and observed values of L q gives new 
impetus for a thermal interpretation of the quiescent emis- 
sion. As we discuss in § spectral fits with realistic emer- 
gent spectra will show an increase in the emitting area over 



one criticism of the thermal emission hypothesis, namely 
that the emitting areas are too small. Of course, other 
luminosity sources, includin g accretion and a shock from 
a rotation-powered pulsar (Campana et al. 1998a), may 
also occur along with thermal emission. An important 
challenge is to identify the luminosity fraction originating 
from the NS surface. 

2. Local Heating in the Deep Crust 

Accretion pushes a given fluid element in the crust 
ever deeper and forces the electron Fermi energy to grow 
wi th time. The electro ns eventually capture onto nuclei 



( Haensel & Zdunik 199C ) , which become progressively neu- 
tron rich and then fuse via pycnonuclear reactions. This 
heat is locally deposited in the crust about the thin layers 
in which the reactions occur. Most of the heat is conducted 
into the core for a steadily accreting cool NS. The interior 
heats up, on a timescale of ~ 10 4 -10 5 yr, until either the 
entire crust is nearly isothermal or, if the accretion rate is 
sufficiently rapid, until neutrino cooling balances the heat - 
ing (Miralda-Escude et al. 1990; [Brown fc Bildsten 1998| ). 

The situation is different for a transiently accreting NS. 
During an outburst, the reactions locally heat the crust 
at a rate QnucM and raise its temperature by ~ 1%. The 
layer in which the reactions occur {shaded region, Fig. |l|) is 
then hotter than both the atmosphere and the core. After 
the rapid accretion stops, the crust cools and its thermal 
profile relaxes to that of an isolated cooling NS, i.e., the 
temperature monotonically increases with depth. While 
the crust heats and cools in response to the changing ac- 
cretion rate, a fraction / of the deposited heat flows inward 
and adds to the heat stored in the core, while a fraction 
(1 — /) flows outward and radiates away. By equating the 
sum of the quiescent luminosity from the cooling core (we 
use the core temperature-luminosity relation for an atmo- 
sphere with an accreted ma ss of 3 x 10" 11 M Q ; jPotekhin] 
Chabrier, fc Yakovlcv 1997|) and the neut rino luminosity 
modified Urea: |Friman fc Maxwell 197E ; crust neutrino 
bremsstrahlung: |Pethick fc Thorsson 1997 ) with the en- 
ergy deposited during the outburst, we determine T c . For 
example, using the core temperature-luminosity relation 
for a fully accreted envelope (Potekhin et al. 1997) implies 
T c « 10 8 (/(M/10- 10 M Q yr- 1 )) a4 K. For these tempera- 
tures, neutrino cooling is unimportant. 

After many outbursts the core reaches an equilibrium 
limit cycle, in which the heat that it loses between out- 
bursts (the incandescent luminosity Li) is replenished by 
the fraction / of the nuclear energy deposited during each 
outburst. We expect that all observed transient NSs will 
have long since reached this equilibrium limit cycle. The 
timescale to heat the core is much shorter than the lifetime 
of a low-mass X-ray binary; at (M) = 10~ n M Q yr _1 , the 
time to heat the core to equilibrium is less than 10 7 yr 
for / ;> 0.01. Since the energy injected into the crust dur- 
ing the outburst period t a reradiates during the recurrence 
interval t r , we expect an average incandescent luminosity 



Li w fQ nuc (i f Mdt 

\*r ./outburst 
= /Qnuc(M). 



(1) 
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Fig. 1. — Timescale for heat to diffuse from a given pres- 
sure to the surface (eq. [gj), for a neutron star (R = 10 km, 
M = 1.4M Q ) accreting at (M) = 10~ n Af Q yr _1 . Most 
of the nuclear energy (98%) is released in the shaded 
region. The change in slope at the shaded region is 
caused by a reduction in the electron to baryon ratio. 
At P > 10 31 ergs cm -3 , the pressure goes from electron- 
dominated to neutron-dominated. 



This relation, although simple, is not useful observation- 
ally because both sides of the equation depend on the 
source distance. To divide out this dependence, we rewrite 
equation (jl]) in terms of fLuences. The outburst fluence is 
L t = t r (M)GM I R and the incandescent fluence is Lit r . 
The ratio of the fiuences is then simply 



L t GM/R 200 



Ljt r 



fQ, 



f 



(2) 



Because the core temperature cannot appreciably change 
in a single outburst, there will be little change in the in- 
candescent luminosity from one outburst cycle to the next. 

While the crust relaxes to a cooling thermal profile and 
radiates the fraction 1 — / of the deposited energy, the 
neutron star luminosity asymp totically approaches Lj o n 
the thermal diffusion timescale ( Henyey fc L'Ecuyer 1969| ) , 



cp 
pK 



1 yr 



1/2 dP 
9 

2 x 10 14 cms~ 2 
9 



P 



10 31 



ergs cm 



3/4 

( 3 ) 



Here P is the pressure, g is the surface gr avity, cp is the 



speci fic he at, and K is the conductivity (see Brown fc Bild , 
sten l |998| for details). In Figure |l| we plot ith as a function 
of pressure for a NS accreting at (M) = 10~ n M yr -1 . 
The thermal diffusion time in the crust is insensitive to the 
temperature and hence is almost independent of accretion 
rate. 

For objects with lower outburst fiuences (e.g., Aql X- 
1, for which t r ~ lyr and t Q sa 30 days), most of the 
heat is stored in the core, i.e., / « 1. The reason is that 
the energy deposited in the outburst only raises the crust 
temperature by - 0.1 (M/10^ 11 M yr" 1 )^,./! yr)%. As 
a result, the luminosity radiated as the crust relaxes is only 
~ 0.3% different from the incandescent luminosity. Note 
that this argument does not hold for objects with similar 



(M) but longer recurrence times, such as Cen X-4 and 
4U 1608-52. For these objects, the crust temperature can 
change by > 10% during an outburst, and so we expect 

/<!• ~ 

Figure g displays L j L q (L q is the observed quiescent lu- 
minosity) as a function of t r /t for some NS (filled squares) 
and black hole (BH) (open squares) transients. Note that 
the BH candidates appear uncorrelated compared to the 
NSs. We also plot the relation for Lj (eq. g) for differ- 
ent values of /. EXO 0748-676 has L a /L q w 1000 and 
from the diagram appears likely to accrete during qui- 
escence. Indeed, it is presently accreting (as indicated 
by the RXTE / AYL-Sky Monitor) at an even higher rate 
(L ~ 10 36 ergss _1 ) than when the L q for Figure || was 
measured. Repeated observations after the NS goes into 
quiescence can discern the fraction of the heat radiated as 
the crust relaxes, i.e., 1 — /. The sum of this fluence and 
the incandescent fluence is just Q nuc {M)t r . These obser- 
vations can in principle constrain g, on which the thermal 
time strongly depends (eq. ||). 
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Fig. 2. — The ratio of outburst luminosity L a to quies- 
cent luminosity L q as a function of the ratio of recurrence 
interval t r to outburst duration t . The lines are for dif- 
ferent fractions, / = 0.1 (dashed line) and f = 1.0 (solid 
line), of Q nuc = 1.0 MeV per baryon mass deposited at 
depths where the thermal time is longer than the outburst 
recurrence time. Also plotted are the observed ratios for 
several NSs (filled squares) and BHs (open squares). For 
most of the BHs, only an upper limit (arrow) to L q is 
known. Data is from Chen et al. (1997), with the e xcep- 
tion of L q for the Rapid Burster (Asai et al. 1996a). For 
Aql X-l and the Rapid Burster, L Q and t D are accurately 
known (RXTE/ All-Sky Monitor public data); for the re- 
maining sources L a and t Q are estimated from the peak 
luminosities and the rise and decay timescales. 
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3.1. Observational Signatures 

The expected heating in the crust naturally generates an 
incandescent luminosity Li that matches the observed val- 
ues of L q . Of course, the star may still a ccrete, either per- 



sisten tly or spasmod ically, in quiescence ( Narayan, Garcia, 
& Mc ^lintock 1997 ). Accretion onto the N S surface at low 



rates will also produce a thermal spectrum ( Zampieri et al , 
19950 so a further challenge is to distinguish between the 
two luminosity sources. If the NS does not accrete during . 

quiesq onco, there will bo (1) little variation in L, from one to that m the puW Be star system PSR 1259-63 flTa- 



of the ion), if detectable with XMM, would allow direct 
measurement of the gravitational redshift. 

Thermal emission does not explain the ha rd power-law 
tail seen in both Con X -4 ( |Asai ct al. 1996b| ) and Aql X-l 
(Campana et al. 1998b). Evidence of an additional lumi- 
nosity source is also indicated by the fast (timescales of 
order days ) variations in the quie scent flux observed from 
Cen X-4 (Campana et al. 1997). If a rotation-powered 
pulsar becomes operational in quiescence, the pulsar wind 
can interact with material in the NS environment, similar 



outburst cycle tn the nc^t (we here mean thr ■increscent vani & Arons 1997). The emission from the shock is most 



luminosity emitted after the crust has relaxed to a cool- 
ing thermal configuration), (2) a slow monotonic decrease 
in flux and T c g while the crust relaxes after the outburst, 
(3) no short-term intensity fluctuations (barring environ- 
mental effects), (4) stability of the ratio of the NS radius 
to distance (as inferred from spectral fitting), and (5) an 
absence of metals in the spectra. 

The quiescent energy spectra of 1608-52 , Cen X-4 (both 
observed with ASCA; [Asai et al. 1996b|), and Aql X-l 



(observed with ROSAT /Posit ion Sensitive Proportional 
Counter; Verbunt et al. 1994 ), when fit by a blackbody 
(BB) with an additional power-law tail for Cen X-4, imply 
substantially smaller radii (1.5 km for 1608-52, 1.8 km for 
Cen X-4, and 1 km for Aql X-l) than expected from a NS 
for the best-fit BB temperatures (0.16-0.30 keV). These 
measured radii have puzzled observers and motivated the 
interpretation that either the quiescent NS luminosity does 
not originate from the surface or only some small fraction 
of the surface radiates. However, the emergent spectra at 
T e ff <J 5 x 10 6 K is very different from a BB. The opacity is 
free-free dominated and is therefore proportional to is~ 3 , 
where v is the photon frequency. As a result, higher en ergy 



photo ns escape from greater de pths, where T > T e g (Fto- 
mani L987; Zampieri et al. 1995). Spectral fits of the Wien 



tail with BB curves then overestimate T c ff and underesti- 
mate the emitting area, by as much as orders of magnitude 
( Rajagopal fe Romani 1996 ; Zavlin et al. 1996). Simple 
comparisons between the observed spectra and the hydro- 
gen atmosphere models imply larger emitting areas in all 
cases. 

An incandescent NS should have a pure H atmosphere 
because the time for h eavy elements to settle out of the 
photosphere is ~ 10s ( Romani 1987 ). However, if accre- 
tion continues in quiescence, there is an accretion rate, Mz, 
above which metals are dumped into the atmosphere fast 
enough to maintain their abundance. In this case a mea- 
surable metal abundance may exist. In the absence of spal- 
lation (which depends on the accretion geometry), metals 
will be underabundant relative to their infalling value for 
M < Mz ss 4 x lO-^AfmVT-^fcRTpfr/ O.l keV) 3 / 2 (8/Z) 
( Bildsten, Salpctcr, fc Wasserman 1992 ). The accretion 



rate Mz coincides with the accretion rate at which the 
resulting luminosity from quiescent accretion, GMM / R, 
is comparable to the thermal emission. // accretion al- 
ters the luminosity, it also alters the photo spheric abun- 
dances. The emergent spectra from an atmosphere with 
solar abundance metal s is detectably different from that of 
a pure H atmosphere (Rajagopal & Romani 1996), mostly 
because of the gravitationally redshifted O vm photoion- 
ization edge. This feature (at 0.87 keV in the rest frame 



likely nonthermal, however, so future observations can dis- 
tinguish the part of the spectrum contributed by thermal 



emission. 



3.2. Aql X-l and SAX J1808.4 3658 

Rece nt observations of Aql X-l as it faded into quics- 
cence (Campana et al. 1998b; Zhang, Yu, & Zhang 1998) 
revealed an abrupt (1 day) decay of the luminosity follow- 
ing a more gradual decline with a roughly 17 day timescale. 
After the sudden dimming, the luminosity persisted at 
L q ss 10 33 ergss^ 1 for the remainder of the observation. 
Both Campana et al. (1998b) and Zhang et al. (1998) in- 
terpreted this sharp transition as the onset of the pro- 
peller effect (centrifugal barrier). Accretion at the low 
rates (« 10~ 13 M Q yr^ 1 ) needed to ex plain L a is very dif - 
ficult if the propeller is operational ( |5tella et al. 1994 ). 
In addition, accretion in spite of the propeller implies a 
flow onto the polar caps that would produce luminosity 
variations at the NS spin period, which are not seen. The 
propeller can be avoided only if the NS in Aql X- l rotate 
with P s > 0.6( J B/10 9 G) 6 / 7 s ( |Verbunt ct al. 1994[ ). Inter- 
preting the 549 Hz oscillation seen during a type I burst 
from Aql X-l (Zhang et al. 1998a) as the rotational fre- 
quency would then imply that B < 10 6 G. This magnetic 
field constraint relaxes substantially if thermal emission 
causes L q and makes more plausible the NS becoming an 
active millisecond radio pulsar. 



The steadiness of L a (over a 20 day observation with 



BeppoSAX , Campana et al. 19981) is naturally explained 
by thermal emission, provided that / « 1 (which we expect 
from the arguments in § 2). The incandescent emission 
from the hot NS has the correct magnitude to explain the 
observations without requiring any ad hoc assumptions. 
Alternate explanations, such as accretion onto the magne- 
topausc or shock emissio n from a rotation-powered pulsar 
(Campana et al. 1998a), are unlikely to produce both a 
steady luminosity and a thermal spectrum, especially one 
with an emitting radius comparable to that of a NS. How- 
ever, at least one other emission mechanism is required to 
account for the hard tail. 

W ith the recent discovery of a 401 Hz accreting pul- 
sar (Chakrabarty fc Morgan 1998 : Wijnands fc van dei 
Klis 1998|) in the transient SAX J1808.4-3658 ( jin't Zanc 
et al. 199q ), an upper bound of 6 x 10 26 G cm s may be 



placed on the dipolc magnetic moment if the magneto- 
spheric radius is presumed to be less than the corotation 
radius. Given a recurrence interval of 1.5 yr, an outburst 
duration of ss 20 days, and an outburst accretion rate of 



3 x lO-^Mgyr" 1 (in't Zand et al. 1998), we expect 



an incandescent luminosity of sa 6 x lO^ergss 1 (//1.0) 
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which corresponds to an unabsorbed flux (4 kpc distance) 
of 3 x 10 -13 ergs cm -2 s -1 (//1.0). The inferred surface 
magnetic field and spin period are sufficient to power 
a rotation-powered pulsar. If this occurs, th e magnc 



tosph cric emis sion is L x < 4 x 10 ergss (Becker & 



Triim jer 1997) and may contaminate the thermal emis- 
sion. 

4. Conclusions 

We have demonstrated that reactions in the inner crust 
of an accreting NS heat the interior enough to make the NS 
incandescent after accretion halts. Transiently accreting 
NSs then radiate ~ 5 x 10 32 -5 x 10 33 ergss -1 in quiescence. 
The quiescent emission is a thermometer that probes the 
heated deep crust an d core. As in the case of isolated 



cooling NSs (see, e.g., [Tsuruta 1998| ), it may be possible 
to test for the presence of enhanced neutrino emissivities 
by placing limits on the interior temperature of the NS. 

Upcoming missions will greatly increase our understand- 
ing of previously discovered transients and open up new 
populations for study. Especially promising arc the low- 
luminosity X-r ay sources in globular clusters (Hertz & 
Grindlay 1983| ). A subject of debate is whether the 



brighter sources of this class (with 10 33 ergss -1 < L x < 
10 34 - 5 ergss -1 ) are quiescent so ft X-ray transients, as op- 
posed to ca taclysmic variables (Vcrbunt, van Paradijs, & 
Elson 1984). If so, future X-ray observations can mine a 
potentially rich source of NS transients. With multiple 
quiescent NSs in a single field of view, all at the same 
distance and reddening, comparisons between sources will 
be easier and radius determinations more certain. Al- 
though in many ways distinct from t he low-mass systems 



the neutron stars in Be trans ients (Apparao 1994; 
den Hcuvel & Rappaport 1987) should also behave as we 
have discussed. These NSs have strong magnetic fields 
(B ;> 10 11 G), show pulsations, and undergo a variety of 
outbursts flBildsten et al. 1997|). 
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